Abstract 



Based on ideas of off-diagonal long range order and two-fluid model, we demon- 
strate that adiabatic phases for a slow motion of a vortex in a superconductor film 
give rise naturally to the Magnus force at finite temperatures. 
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There is an explosive renewed interest in vortex dynamics in superconductor sflX|. A 
consistent and clear understanding of the basic physics for vortex dynamics is needed. It 
has been shown recently by the aid of the many-body wavefunction at zero temperature 
that adiabatic phases for a slow motion of vortex line in a superconductor give rise natu- 
rally to the Magnus force 0. Here we demonstrate that the same result can be extended 
to finite temperatures, with the many-body wavefunction replaced by the density matrix. 
For a clear representation, we consider a homogeneous 2-d superconductor film taken to 
be the x — y plane. As in Ref . , we first discuss the case that the magnetic field effect 
is negligible. 

The density matrix for N electrons is p^ = PN( r i--- r N', r 'i--- r jv) with {r?} the position 
vectors and spins of the N electrons. The one particle reduced density matrix is defined 

as 

N 

pi(r; r') = N Y[ d 2 r j p N (r...r N ; r', r 2 ...r N ). (1) 

J j=2 

For a superconducting state, a two-fluid model may be applied for the one particle density 
matrix The essence of the two-fluid model is that, at any instant there exist two 

real functions, p„(r, r') and p s (r, r') such that[|j] 

pi(r; r') = p n (r; r') exp{z[r • k n (r) - r' • k n (r')]} 

+ Ps (r;r')exp{*[0(r)-0(r')]}. (2) 
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Here real functions p n , p s , k n , and are the functions of time and temperature. The total 
electron number density is p(r) = p n (r;r) + p s (r;r). If p s = 0, this would describe the 
system with classical velocity field v n = Kk. n /m with m the effective mass of an electron, 
and the system is in the normal state, where p n — > exponentially as |r — r'| — > oo. If 
p n = 0, the system is at zero temperature as discussed in Ref.0. We have assumed that 
the wavevectors k n and V r 0(r) do not alter appreciably over a distance of the order of the 
coherence length £. The existence of finite superfluid density p s at finite temperatures is 
the manifestation of off-diagonal long range order, but we cannot write p s as the product 
of wavefunction[[|]. Instead, it should be determined from the pair wavefunction in two 
particle reduced density matrix because electrons form pairs. 
The Berry phase for a pure state is@, 

e = ifdr - < * B (r )|Vro*«(ro) > . (3) 
Writing this equation in a more explicit and suggestive form, the Berry phase becomes 

N 

6 = i j> dr • Jim V ro / Y[ 

Jr r o — >r o J j—i 

PN,v({ r j}> r o;{ r j},*'a)- ( 4 ) 
Here p^ tV is the density matrix defined from the many-body wavefunction for a vortex ty v 
as pn,v = ^({r}, ro)^^! 1 "'}, r o)- The vortex position is r . We do not know the exact 
form of the Hamiltonian for vortex dynamics. Instead, we have a good knowledge of the 
wavefunction (or the density matrix) describing a vortex state, which is the starting point 
to derive the Magnus force by calculating the adiabatic phases. This method is essentially 
the Born-Oppenheimer approximation where the vortex coordinates are slow ones, and 
has been put into a more general perspective with the aid of the Berry phase ]/]]. Eq.(4) 
is the generalization of the Berry phase from a pure state to a mixed state. 

In the region far away from the vortex core where the supercurrent generated by the 
vortex is small, the two-fluid model can be applied. The ansatz for the one particle reduced 
density matrix of the vortex state at a finite temperature can be written as eq.(2), with 

TYl 

0(r) = 0(r - r ) + — v s ■ r , (5) 

where 9 = q v arctan[(y — yo)/(x — xq)]/2 with q v the winding number of a vortex and v s 
is the superfluid velocity. 

Using eqs.(4), (1), (2), and (5), the Berry phase for a vortex adiabatically moving 
around an arbitrary close loop V is 

6 = - jf rfr / d 2 rV ro £(r - r Q )^p) , (6) 
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where there is no contribution from the normal fluid p n . Ignoring the contribution from 
the core as done in Ref.[@, we obtain the Berry phase 

6 = -2-Kq v ^S(V) . (7) 

The dynamical phase corresponds to the potential energy of the system at a finite 
temperature. For given density matrix, the potential energy of the system is given by 
U(r ) = J JJjLi dvj lim r /^ rj . H({rj})p N:V ({vj}, r ; r'-, r ). One and two particle reduced 
matrices are involved for the usual BCS Hamiltonian H. The one particle reduced density 
matrix is given by eq.(2). For the two particle reduced density matrix it may be reasonable 
to assume that this density matrix is close to the one without the vortex, except near the 
vortex core where the contribution can be put into the core energy. Then the potential 
energy depending on the position of the vortex in the presence of an external current is 

U(t ) = -^-hq v r ■ v s x z . (8) 

There is no dependence on the normal current. 

Combining eqs.(7) and (8), we obtain the Magnus force at a finite temperature: 

F m = y%,(v s - r Q ) x z , (9) 

which is the same as eq.(9) in Ref. @, except that here the superfluid electron number 
density is temperature dependent. 

As demonstrated in Ref. 0, the electromagnetic field associated with a vortex does 
not affect the Magnus force because of the charge neutrality of the superconducting film. 
The disorder does not affect its existence either, except for a possible reduction of the 
superfluid density due to the localization of some Cooper pairs and a pinning force due 
to the introduced inhomogeneity. 

Finally, we wish to point out that neither the Iordanskii force || nor the conservative 
force recently discussed by Ferrell[|] seem to appear in the present approach. 

We thank M.Y. Choi and Y. Tan for useful discussions. This work was supported by 
US National Science Foundation under Grant No. DMR 89-16052. 
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